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Abstract. We write the Mathisson-Papapetrou equations of motion for a 
spinning particle in a stationary spacetime using the quasi-Maxwell formalism 
and give an interpretation of the coupling between spin and curvature. The 
formalism is then used to compute equilibrium positions for spinning particles 
in the NUT spacetime. 



Introduction 

The purpose of this paper is to use the quasi-Maxwell form of Einstein's held 
equations for a stationary spacetime to give an interpretation of the force that acts 
on a spinning particle according to the Mathisson-Papapetrou equations. 

The paper is divided into hve sections. In the hrst section, we briefly review the 
quasi-Maxwell formalism for stationary spacetimes. The Mathisson-Papapetrou 
equations of motion for a spinning particle, as well as the various spin supplemen- 
tary conditions, are discussed in the second section. The third section recalls the 
main formulae for a magnetic dipole on a magnetic field, which are compared with 
their gravitational counterparts in the fourth section. It turns out that there is a 
remarkable analogy between the two sets of formulae, shedding light into the na- 
ture of the spin-curvature coupling. In the fifth section, as an example, we use the 
quasi-Maxwell formalism to compute equilibrium positions for spinning particles in 
the NUT spacetime. 

We use the notation and sign conventions in |Wal84] . where latin indices repre- 
sent abstract indices, except for the latin indices , which will be numerical 
indices referring to the space manifold (i.e. = 1, 2, 3). As is usual, we will 
not worry about the vertical position of space indices on orthonormal frames. 

1. Quasi-Maxwell formalism for stationary spacetimes 

We start by briefly reviewing the quasi-Maxwell formalism for stationary space- 
times. For more details, see fOhOl ICN05I IEmb841 ILL971 ILBNZ98I INZ971 1537)0] . 

Let (M, g) be a chronological stationary spacetime with a complete timelike 
Killing vector field T. Then the action of R on M by the flow of T makes it a trivial 
principal R-bundlc [AndOO, BSOO , and the quotient space M/R can be identified 
with a 3-dimensional submanifold E C M. One can think of this quotient as the set 
of all stationary observers. If 7T : M S is the quotient map, p £ £ and v € T p T,, 
we define to be the vector field along the integral line of T through p which 
is orthogonal to T and such that tt*v* = v. We can then define the Riemannian 
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metric h on £ through 

h(v, v) — g{v\ it). 

This metric has the physical interpretation of yielding the distance measured be- 
tween nearby stationary observers through Einstein's light signaling procedure 
|LL97j . The space manifold is the Riemannian manifold (£,/i). 

There exists a function <f> : £ — > R, called the gravitational potential, such 
that 

g{T,T) = - e 2 W 0lr ). 
The gravitational field is then the vector field G defined on £ by 

G = grad (f>. 

It is easy to show that — G^ is simply the acceleration of the stationary observers. 

For a particular choice of the submanifold £ C M, one can define the one-form 
A through 

A(v) = -e-^g(T,v) 

for all v £ T£ C TM . It turns out that although A depends on the choice of S, 
the 2-form 

H = -e^dA 

does not, and hence is well defined on the quotient. Assuming that £ is orientable, 
we define the gravitomagnetic field to be the vector field 

H = (*H) S , 

where * is the Hodge star and " is the usual map between vectors and covectors 
determined by h. It is not hard to show that jH7 is simply the vorticity of the 
stationary observers. 

The names above are justified by the fact that the projection of any timelike 
geodesic on the space manifold, parameterized by the proper time r, satisfies the 
equation 

— - = 7 2 G + 7 u x H, 

ar 

i 

where u is the tangent vector and 7 = (l + u 2 ) 2 is the energy per unit mass of 
the corresponding free-falling particle as measured by stationary observers. As a 
consequence of this equation, the particle's total energy per unit mass 

E = e 7 

is conserved. 

If {Ei, E2, E3} is a right-handed orthonormal local frame on the space manifold, 
then {E , E\,E\, e\} is an orthonormal local frame on M, where Eq — e~^T. The 
nonvanishing Christoffel symbols associated to this frame are Oli02 
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where Tj k are the Christoffel symbols associated to {E%, E2, E3}. The components 
of the Ricmann curvature tensor are therefore given by 



^RoMj — — VjGi + GiGj — -HikHkj; 
^Roijk = - (VjHik — VkHij — 2GiHjk) : 



2 

^Rijki — Rijki + {^HijHki + HikHji — HaHjk) , 
and the components of the Ricci curvature tensor are 
(4) i? 00 = -divG + G 2 + iH 2 ; 
WE i = ~(curlH).-(GxH) i; 

WjZ« = Rij + ViGj - G,G, + - 2 11,11, + ^H"//,,. 
In particular, one can write Einstein's equations in the quasi-Maxwell form 

divG = G 2 + ^H 2 - 4tt (p + tr<r) ; 
curlH = 2G x H 167rj; 

Ric + VG» = ® G» + ® H" — ^H 2 fr + 87tct + 4tt (p - tr ct) fc, 

where Ric is the Ricci curvature of the space manifold and p — T 00 , j = T 0z Ei and 
a = TijEf <X> Ej are the energy density, the energy current and the spatial stress 
tensor measured by the stationary observers. These equations imply that the scalar 
curvature of the space manifold is 

R= 167T/5- ^H 2 . 

2. Spinning particle in General Relativity 

The Mathisson-Papapetrou equations [Mat371 |Pap51| for the motion of a spin- 
ning particle, obtained by neglecting all multipoles higher than mass monopole and 
spin dipole lDix701 EDSU3] . are 

VjjP a + ^R a bcd U b S cd = 
VuS - P AU = 

where U is the unit tangent vector to the particle's history, P is its energy- 
momentum vector (not necessarily parallel to U) and S is a rank-2 antisymmet- 
ric tensor representing the particle's angular momentum. This system of equa- 
tions is underdetermined, and must be closed by means of a spin supplemen- 
tary condition. The usual choice, justified in part by the uniqueness results in 
|Bei671 ISch79al ISch79bj . is the Tulczyjew-Dixon condition |Tul59l 10570] 

S ab P b = 0. 

Less popular choices include the Mathisson-Pirani condition [Mat37[ Pir56 
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(for recent reviews of these matters see for instance |dFC95[ ISem99j ; recent work 
on solutions of these equations in Schwarzschild and Kerr spacetimes can be found 
in [SM961 ISirH99l IHar031 lBdFG04bl lBdFG04al ISm05l IMS06] ). 

Dixon's supplementary condition implies that the angular momentum tensor is 
of the form 

nab \ abed n 



r a p6 ncd 
bcd r ° 



where e is the volume element and m = y/—g(P, P) is the particle's dynamical 
mass (which can be seen to be constant). The spin vector 

= " ^ 
2m 

(which is the relativistic analogue of the angular momentum vector) is orthogonal 
to P (hence spacelike), and satisfies 

m z 

It is possible to show [K72l ITdFC761 ISem99j that 

(1) U a = (pa + }_ S *b Rbcde pc s de 

where 

A = 2m 2 + X -RabcdS ah S cd . 
Therefore U and — P differ bv a Lorentz transformation of order 



where r is the local radius of curvature. Assuming that the particle is not ro- 
tating too fast, we can take P and U to be parallel, P — mil, in which case 
Dixon's and Mathisson's spin supplementary conditions coincide. The (approxi- 
mate) Mathisson-Papapetrou equations are then simply written as 

(2) mwer - A R - bcdU b e ^f UeSf . 

(3) V V 8=g{V v U,a)U. 

These equations are consistent with the spin supplementary condition g(U, s) = 0, 
as s is Fermi- Walker transported along the motion. 

The first equation indicates that a net force arises from the particle's spin cou- 
pling to the spacetime curvature. The expression for this force, however, is not 
particularly enlightening. We will see that within the quasi-Maxwell formalism its 
expression is actually quite natural. 

3. Magnetic dipole on a magnetic field 

As is well known (see for instance |Jac98] ). a magnetic dipole with moment m 
placed in a magnetic field B will experience a torque 

N = m x B, 

which can be thought of as arising from the potential energy 

U = -m- B. 
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Such dipole will experience a net force 

F = VB m (div B) m, 
where the last term is not usually written because of Maxwell's equation div B = 0. 

4. Spinning particle in Quasi-Maxwell formalism 

We will consider the special case when the particle is at rest with respect to the 
stationary observers, i.e. 

U = E . 

The spin supplementary condition then yields 

s° = 0. 

Using the formulae for the Christoffcl symbols, it is easily seen that the equation 
for the Fermi- Walker transport of s can be written as 

^f! + (4) T i «j _ o ^ — - iff s J - 
dr + 0J dr 2 v ~ ' 

i.e. 

Da 1 

^ 1F = - 2 Hxs - 

This formula shows that the particle's angular momentum changes under a gravito- 
magnetic field exactly like a magnetic dipole moment under a magnetic field (except 
for a factor ~). 

Using the expressions for the components of the Riemann tensor, one can com- 
pute the force due to coupling between spin and curvature to have components 

F i = ~\ {4)Rl ojk^i-l)*! = \ (Vi^ife - Vfcffy - 2GiH jh ) e jH s l 
where e is the volume element of the space manifolcQ. Therefore 

F % = — (VjHik ~ GiHjk) Sjkis 1 = — (V ' j{eik m H m ) — Gi£jk m H m ) Ejkis 1 



- \ (5 l3 8 ml - 6 u S mj ) (VjH^s 1 - l - (S 3] S ml - 5 3l 5 m3 ) G^s 1 
= \{ViH l )s l - \{V 3 W)s l - l - (3£ roi - 8 ml ) G^s 1 
= \{ViHl)j - \{^ 3 Wy - (ffVjGi, 

that is, 

(5) F = ivH ■ s- - (divH)s- (s ■ H) G. 

Therefore the spin-curvature coupling can be interpreted as consisting of a grav- 
itomagnetic part, which mimics the magnetic force on a magnetic dipole (apart 
from the factor |), plus the weight of the dipole energy. 

This similarity between the spin-curvature coupling and the force on a magnetic 
dipole was noticed long ago in the context fo the linearized theory by Wald |Wal72) . 
Equation ([5]) extends this analogy to an arbitrarily strong (albeit stationary) field. 
Notice that only the first term in ([5]) appears in |Wal72j . as the other two are not 
linear in the fields (it is easily seen that divH = — G • H). 



^We have e ^* = —£ijk on anv orthonormal frame. 
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5. Spinning particle in NUT spacetime 

We would like to use equation ([5]) to compute equilibrium positions for spinning 
particles in stationary spacetimes. Unfortunately, this is a situation where the 
approximate equations (J2]) and (|3|) are expected to break down. Indeed, since the 
components of the Riemann tensor are quadratic in the components of the fields G 
and H, we expect these to be of order -, where r is the local radius of curvature. 
Therefore the force on the spinning particle is of order 

l 4- 

For the particle to remain at equilibrium, this force must balance the gravitational 
force mG, which is of order — . Therefore, we expect 

N ™> „ \s\ 1 

r z r TJIT 

at any equilibrium position, whereas the approximate equations assume that this 
quantity is small. 

Equation (fTJ), however, implies that the relation U = is exact if s x F = 0, 
where F is the force on the spinning particle given in ([5]); for an equilibrium position, 
this is equivalent to requiring that s x G = O0 Therefore the equilibrium positions 
computed using equation ([5]) will be exact for the Tulczyjew-Dixon supplementary 
condition if the angular momentum is aligned with the gravitational field. 

For completeness, let us remark that the Mathisson-Pirani spin supplementary 
condition leads to the motion equation 

1 
2 

where the constant m is now defined as m = -g(P, U) |Mat37| IPir56| IMTW73| 
ISte04j . For an equilibrium position, this will reduce to the approximate equation 
<j2j) when s x (G x H) = 0. Therefore the equilibrium positions computed using 
equation ([5]) will be exact for the Mathisson-Pirani supplementary condition if the 
angular momentum is aligned with the "Poynting vector" G x H. Thus for instance 
the equilibrium positions along the axis of the Kerr-de Sitter solution, computed 
in |SK06| using the Mathisson-Pirani supplementary condition, coincide with the 
equilibrium positions computed using Tulczyjew-Dixon supplementary condition, 
and can both be obtained from ([5]). 

As an example, we now compute equilibrium positions for spinning particles in 
NUT spacetime. Recall that the NUT spacetime is a stationary solution of the 
vacuum Einstein field equations describing a gravitomagnetic monopole, given in 
local coordinates by 

g = -e 2,p (dt + 2lcos9dtp) 2 + e^dr 2 + (r 2 + l 2 ){d6 2 + sin 2 Odtp 2 ), 

where 

Mr + I 2 



mVuU a = -^R a bcd U b S cd + S ab VuVuU b , 



e 



20 _ 1 



r 2 +P 

and M, I are two parameters representing the total mass and (half) the gravito- 
magnetic charge |NTU63| ILBNZ98} ISKM03] . The metric of the space manifold 



2 It is interesting to note that the linear momentum of the field produced by a magnetic dipole 
m on an electric field E has the approximate expression Pfldd = E X m Jac98 . 
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determined by the stationary observers is 

h = e^dr 2 + (r 2 + l 2 )(d0 2 + sin 2 9dtp 2 ), 
and the gravitomagnetic potential 1-form is 

A = 21 cos 9dip. 
Thus the gravitomagnetic held corresponds to the 2-form 

H = -e^dA = 21^ sin OdO A dip, 

and hence to the 1-form 

H * =* H = -TTn dr > 

implying that H is radial. If we choose the angular momentum s to be also radial, 
then the equilibrium positions will be exact for the Tulczyjew-Dixon supplementary 
condition. Moreover, since the 1-form associated to the gravitational field is 

G i = -d<p = -4>'dr, 

where the prime denotes differentiation with respect to r, we have G x H = 0, and 
the equilibrium positions will also be exact for the Mathisson-Pirani supplementary 
condition. 

By equation (j4]), the angular momentum s will be constant. We now compute 
the 1-forms associated to each term on the expression ([5]) for the force on a spinning 
particle (here we must keep track of the vertical position of the indices, as we are 
working with a non-orthonormal coordinate basis). To begin with, we have 

(VH • s) J = V r H r s r dr = (H' r - T r rr H r )s r dr = {H' r + cf>'H r )s r dr 

(the relevant Christoffel symbols can be readily computed from the geodesic equa- 
tions of the space manifold). Moreover, 

(divH)s B = -(G ■ H)s r dr = -G r H r s r dr = -G r H r s r dr = <p'H r s r dr. 

Finally, 

(s • H)G J = -s r H r <P'dr. 
Collecting these terms together yields 

F* = ^(H' r + 2cj) , H r )s r dr. 
The equation for the equilibrium positions is then 

s r (e 2 ^)' 



mG« + F» = <=> -2m<j>' + (H' r + 2cj) , H r )s r = <=> 



m {H r e 2r t>) r 



For r 3> v M 2 + I 2 , we have 



and hence the equilibrium condition reduces to 

s r r Mr 
m ~ W 

In particular, the angular momentum must point towards the center if I > 0. This 
could easily be computed from the linearized formulae in jWa!72] . In the strong 
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field region, however, things are more complicated. For instance, near the horizon 
r H = M + \/M 2 + I 2 one has 

and hence the equilibrium condition reduces to 

m ~ H r {e 2 *)' ~JT r ~ 21 ' 

In particular, the angular momentum must point away from the center in this 
region. 
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